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ABSTRACT 

The morphology of the e Eridani dust ring is reproduced by a numerical simulation of dust particles 
captured into the 5:3 and 3:2 exterior mean-motion resonances with a 0.3 eccentricity 10 -4 solar mass 
planet at periastron at a semi-major axis of 40 AU. The morphology will differ when the planet is 
at apastron, in about 140 years. Moderate eccentricity planets in outer extra-solar systems will cause 
observable variations in the morphology of associated dusty rings. 

Subject headings: Subject headings: stars: circumstellar matter — stars: planetary systems — stars: 
individual (Epsilon Eridani) 



1. INTRODUCTION 

The Infrared Astronomical Satellite (IRAS) measured 
far infrared fluxes from the vicinity of main sequence stars 
(Aumann et al. 1984), including the nearby systems Vega 
and Formalhaut which are of moderate age, 10 8 years, at 
the end of the era when rocky planets are expected to 
form. Recently, submillimeter wavelength images of stars 
with infrared excesses have revealed structure in some of 
these disks. Formalhaut has a dust ring with a central cav- 
ity (Holland ct al. 1998), Vega, (3 Pictorus and HR 4796A 
have secondary dust emission peaks and central cavities 
(Wilner et al. 2002; Holland et al. 1998; Jayawardhana et 
al. 1998; Koerner et al 1998), HR 4796A has a ring ev- 
ident in scattered light from the central star (Schneider 
et al. 1999), and e Eridani has a 60 AU dusty ring with 
peculiar 4 peaks of emission (Greaves et al. 1998). HR 
4796A and j3 Pictorus are younger, at ages of only ~ 10 7 
years, whereas e Eridani (HR 1084) is a nearby 0.5-1 Gyr 
old K2V (0.8M Q ) star (Soderblom & Dappen 1989; Song 
et al. 2000). 

Because of the short lifetime of the dust particles in 
these disks, and the scarcity of gas, the dust must continu- 
ally be replenished from collisions of larger orbiting bodies, 
possibly analogous to Kuiper belt objects. Hence the disks 
in the older systems are commonly called debris disks. 
Dust particles in our solar system spiral toward the sun due 
to radiation drag forces, known as Poynting-Robertson (P- 
R) drag and interactions with the solar wind (Burns et al. 
1979). Gravitational perturbations from planets can scat- 
ter the particles or capture them into orbital resonances 
with planets (e.g., Dermott et al. 1994). Since the life- 
time of captured particles can greatly exceed the lifetime 
of particles which are not captured, resonant trapping can 
sculpt the dusty disk. 

Structure predicted in the dust distribution in the 
Kuiper Belt is primarily a result of resonance trapping 
(Liou & Zook 1999; Moro-Martin & Malhotra 2002). 
Wilner et al. (2002) suggested that the two peaks seen at 
1.3mm in the Vega disk could be explained by the capture 
of dust particles into the 2:1 exterior mean-motion res- 
onance with a high eccentricity three Jupiter mass giant 
planet at a semi-major axis of 40AU. Ozernoy et al. (2000) 



proposed that the morphology of the submillimeter images 
of the e Eridani system could be caused by dust particles 
trapped in high libration orbits associated with the 3:2 
mean motion resonances with a ~ 0.2 Jupiter mass planet 
in a circular orbit. While Ozernoy et al. (2000) considered 
the morphology caused by dust particles carefully placed 
in individual resonances, here we initially distribute dust 
particles over a range of semi-major axes and eccentric- 
ities and search for stable dust configurations which can 
account for the four peaks observed in the submillimeter 
images of the e Eridani system. 

2. NUMERICAL INTEGRATION 

We numerically integrate the orbits of dust particles us- 
ing a conventional Burlisch-Stoer numerical scheme. The 
forces on the particles are gravity from a central star and 
one planet and P-R drag from the central star. The motion 
of the dust particles was restricted to the plane containing 
the planet and star. The ratio of the force from radiation 
to that from the gravity of the central star is parametrized 
with the dimensionless parameter, /?, which depends pri- 
marily on particle size (Burns et al. 1979). Length scales 
are given in units of the perturbing planet semi-major axis, 
a p , and timescales in units such that GM, = 1 where iW* 
is the mass of the central star. The ratio of the planet 
mass to the stellar mass is denoted as //. 

Dust particles see a reduced radial force from the cen- 
tral star because radiation pressure opposes gravity. This 
causes the location of the mean motion resonances with the 
planet to shift radially by a factor of (1 — /3) 1 / 3 . However 
the magnitude of the P-R drag force caused by relativistic 
effects depends on the ratio of the particle velocity to the 
speed of light, c, times (3. Because of our choice of units, 
the speed of light, 



adopted c/ y/GM*/c 



is in units of y/GM*/a p . We have 
= 5 x 10 4 which corresponds to hav- 



ing a stellar mass of 1M Q and a planet radius of 25AU. 
Because of the weak dependence of c on stellar mass and 
planetary semi-major axis, our simulations should scale to 
outer solar systems, such as the e Eridani disk, with the 
caveat that the particle drift rates should be a factor of a 
few times lower. 

With the exception of Wilner et al. (2002), most simu- 
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lations of dust dynamics in planetary systems have been 
carried out for low eccentricity systems. As Wilner et al. 
(2002) found, high eccentricity planets can cause strong 
azimuthal density contrasts in the distributions of objects 
trapped in outer resonances such as the 2:1 and 3:1 mean 
motion resonances. Typically we find that the azimuthal 
contrasts are lower in the outer resonances than those 
closer to the planet. Previous simulations have predicted 
that resonances close to the planet such as the 4:3, 3:2, 
5:3, and 7:4 resonances efficiently trap and hold dust par- 
ticles in Kuiper belt (Liou & Zook f999; Moro-Martin & 
Malhotra 2002). We restrict our simulations to the region 
containing these resonances. 

We began the simulations with 300 dust particles dis- 
tributed randomly in semi-major axis between 1.1 and 1.5 
times the planet's semi-major axis. Longitudes of peri- 
astron and mean longitudes were chosen randomly. Dust 
particle eccentricities were chosen to lie between zero and 
0.4. We ran each simulation for 2 x 10 4 planetary orbit 
periods. The positions of all the particles were recorded 
every planetary year, consequently the phase of the planet 
was the same during each output of the integration. We 
binned the particle positions counting the positions the 
particles at each output of the simulation. The resulting 
spatial distributions and semi-major axis distributions are 
shown in Figure 1 for the planet at periastron and at apas- 
tron. All dust particles contributed to the time integrated 
dust distributions, however, particles which were quickly 
ejected from the system, or quickly impacted the star or 
planet, contributed less to the integrated dust distribu- 
tions. Particles which remained trapped in resonances for 
long periods of time contributed substantially to the dust 
distribution. 

In Figure Id we show the particle distribution restricted 
to semi-major axis between 1.2-1.3, containing the 3:2 
resonance, and between 1.3-1.37, containing the 5:3 res- 
onance, for the simulations shown in Figure la,b when the 
planet is at periastron and at apastron. Only one family 
of resonant orbits dominates for both resonances and the 
phase of the density peaks with respect to the planet de- 
pends on the planet's orbital phase. The distribution of 
particles in the 3:2 resonance exhibits a three-fold sym- 
metry and that of the 5:3 exhibits a five-fold symmetry. 
When the planet is at periastron, the sum of the particle 
distributions in the 3:2 and 5:3 resonances exhibits four 
peaks similar to those observed in the e Eridani system. 

To simulate the effect of the decreasing radiation field 
the dust distribution was weighted by r~ 2 where r is its 
radial distance from the star. The resulting predicted in- 
tensity image for a planet at periastron, and smoothed to 
the approximate resolution of the submillimeter observa- 
tions is shown in Figure 2 with the submillimeter image 
by Greaves et al. (1998) for comparison. The dusty disk 
is tilted assuming a disk inclination of 40° (where 0° cor- 
responds to face on). This is somewhat higher than the 
25° estimated by Greaves et al. (1998) from the ring but 
consistent with the 30° ± 15° estimated from optical line 
data Saar & Osten (1997). 

The morphology of the dust simulation shown in Fig- 
ure 2 exhibits two dominant peaks at mean longitudes of 
approximately ±37r/5 from the planet which correspond 
fairly well to the two dominant peaks south east and west 



of the star observed in the submillimeter map. Near the 
planet, a weaker arc is displayed in the simulation, which 
corresponds to the weaker arc of emission observed to the 
north of the star in the submillimeter map. However, the 
dust emission predicted opposite to the planet is more dis- 
tant from the star than that observed south of the star in 
the submillimeter image. 

In Figure Id we show particle distributions for restricted 
ranges of semi-major axes. When the planet is at perias- 
tron, density enhancements caused by the 3:2 resonance 
add to those from the 5:3 causing the two dominant peaks 
seen in Figure la. When the planet is at apastron, the con- 
centrations caused by the 3:2 resonance are out of phase 
with those in the 5:3. Some substructure is also caused by 
the 7:4 resonance in Figures la,b. 

In Figures 1 and 2 we have given axes in units of the 
planets semi-major axis. The radius of the ring was esti- 
mated to be approximately 60 AU which would correspond 
to a planet semi- major axis of about 2/3 this or 40 AU. 
Assuming a stellar mass of 0.8 Mq, the planet orbital pe- 
riod is 280 years. The total time of our simulations (2 x 10 4 
times the planet period) corresponds to 6 x 10 7 years. 

Dust particles can be trapped in a k + q : k exterior prin- 
ciple mean-motion resonance with resonant arguments of 
the form 

<j> = (k + q)X - k\ p - pw + (p— q)w p (1) 

where p is an integer between and q, A and X p are the 
mean longitudes of the particle and planet respectively, 
and w and zu p are the longitudes of periastron. We con- 
sider the particle to be trapped in the resonance when <fi 
librates about a particular fixed value. Take the planet's 
longitude of periastron as a reference angle (zu p = 0). In 
resonance the resonance angle <f> « mr where n is odd if q 
is even and n is even otherwise. 

For the mean longitude of a dust particle very close to 
the planet to be strongly locked to that of the planet, the 
resonance responsible is likely to be the p = one. In this 
case 



For the 3:2 resonance and when the planet is at periastron 
(X p — uj p = 0), we find dust particles near A = 0, ±27r/3. 
When the planet is at apastron (A p = ir), A is the same. 
This only looks out of phase in Figure Id because we have 
rotated the images so that the longitude of the planet is 
fixed. When the resonant argument <f> is fixed, we will only 
find particles with these longitudes, causing the concentra- 
tions we see in Figure 1. For the 5:3 resonance, when the 
planet is at periastron, A = n, ±37r/5, ±7r/5, and when the 
planet is at apastron, A = 0, ±27r/5, ±47r/5. This five-fold 
symmetry is also consistent with the concentrations seen 
in Figure Id. 

In our simulations we also found particles populating 
resonances with other values of p, however these tended 
to have shorter lifetimes and so contributed less to the fi- 
nal particle distributions. When p ^ 0, higher particle 
densities arise at angles when the particles reach apastron 
or periastron (see discussion in Wilner et al. 2002). For 
example, p = 1 is favored for the 3:2 resonance when the 
planet is at low eccentricity causing two peaks in the par- 
ticle distribution rather than three, which is what is seen 
in the simulations of (Liou & Zook 1999). 
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3. SUMMARY AND DISCUSSION 

We have found that the morphology of dusty mate- 
rial trapped in exterior resonances with a planet can be 
strongly dependent upon the eccentricity and orbital phase 
of the planet. The morphology of the e Eridani dust ring 
seems to be reproduced by dust particles captured into the 
5:3 and 3:2 exterior resonances with a moderate eccentric- 
ity e p ~ 0.3, planet near periastron. When the planet 
was more massive (/j>5x 10~ 4 ) or more highly eccentric 
(greater than 0.4), the resonances closest to the planet did 
not capture and hold particles as efficiently. When the 
planet mass was below a Saturn mass, the 4:3 resonance 
also contributed to the particle distribution resulting in 
dominant symmetric four peaks. The asymmetries ob- 
served in the e Eridani image were not observed in the 
dust distribution so we consider planets less massive than 
Saturn unlikely to explain the particle distribution. When 
the planet eccentricity was below 0.15, the azimuthal den- 
sity variations in the dust distribution were too low to 
account for the morphology in the e Eridani disk. 

Our model differs from that proposed by Ozernoy et al. 
(2000) in a number of ways. The period of the orbital 
planet in our model is 280 years which should cause the 
pattern to revolve by about the star by ~ 1.3°yr _1 , faster 
than the ~ 0.7°yr _1 estimated by Ozernoy et al. (2000), 
where the model planet semi-major axis is ~ 60AU rather 
than at ~ 40AU. Our model planet is located to the north 
of the star, rather than to the west of the star. Our model 
planet mass is similar to theirs but at a moderate eccentric- 
ity. Our model dust concentrations are a result of segrega- 



tion in the phase of a resonant angle, rather than caused 
by a large libration amplitude. Furthermore, because of 
the eccentricity of the planet, our model predicts that the 
morphology of the dusty ring will vary, as well as revolve 
as the planet orbits about the star. 

The initial conditions of our simulations cause many of 
the particles to begin trapped in the resonances and at 
fairly low eccentricity. It is possible that low eccentricity 
planctesimals exist in the system and that they are the 
source of the dust particles that we see. Alternatively the 
planctesimals in the system are further out and the dust 
particles become trapped in the resonances as they spi- 
ral in toward the planet. Further simulations would be 
required to differentiate between these possibilities. In ci- 
ther case, resonant capture into these resonances is less 
likely and less prolonged when the planet mass or eccen- 
tricity is high. If the planet mass is too low then it cannot 
be responsible for clearing a gap or central region in the 
dust distribution. Additional and more massive planets 
would be required to do this. 

It is possible that high eccentricity planets are common 
in the outskirts of extra-solar systems. If so then the re- 
sulting dust distributions would not only revolve (Ozernoy 
et al. 2000), but will also be dependent upon the orbital 
phase of the planets. This is an exciting prospect because 
there would be variations in the dust morphology on ob- 
servable timescales. 

This work would not have been carried out without help- 
ful discussions with Stephen Thorndike, Joel Green and 
Dan Watson. 
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Fig. 1. — a) The dust distribution for a planet at periastron with ratio of planet mass to stellar mass = 10~ 4 , and eccentricity of e p = 0.3. 
The dust particles have /3 = 0.1. The star and planet are denoted as black dots at the origin and at a radius of 0.7 from the origin, respectively. 
Axes are given in units of the planet semi-major axis, b) Same as a) but the planet is at apastron. c) Semi-major axis distribution for the 
particle distribution shown in a. d) Particle distributions are shown for limited ranges of semi-major axes for the simulations displayed in 
Figures la,b. Particles in the 3:2 resonance are shown for the planet at periastron (top left) and at apastron (bottom left). Particles in the 
5:3 resonance are shown for the planet at periastron (top right) and at apastron (bottom right). 
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Fig. 2. — On the right is shown the simulated intensity distribution based on the simulation shown in Figure la, for a 0.3 eccentricity planet 
at periastron of planet to stellar mass ratio 10~ 4 . For comparison we show on the left the 850/im submillimctcr emission map of the e Eridani 
system, (this is Figure 1 by Greaves et al. 1998). The location of the star is denoted in both panels as a star. For the simulated intensity 
distribution we have chosen a disk inclination of 40° where 0° is face on. The dust distribution has been weighted by r~ 2 to simulate the 
intensity field from the star and then smoothed to the approximate resolution of the submillimeter images. We estimate that the planet is 
currently located to the north of the star (denoted as a filled circle in the right hand panel) at a declination of about — 9 d 27 m 20 s (epoch 1998 
as from Figure 1 of Greaves et al. 1998). Axes for the simulation are given in units of the semi- major axis of the planet. 
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